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Abstract 



Recently, it has been possible to construct single-electron transistors to 
study electronic properties, including superconductivity, in metallic grains 
of nanometer size. Among several theoretical results are suppression of su- 
perconductivity with decreasing grain size and parity effect, i.e. dependence 
on the parity of the number of electrons on the grain. We study how these 
results are affected by degeneracy of energy levels. In addition to the time- 
reversal symmetry, for certain energy spectra and more generally for lattice 
symmetries, energy levels are strongly degenerate near the Fermi energy. For 
a parabolic dispersion, degeneracy d is of the order of kpL while the typical 
distance between the levels is of the order of ep/(kpL) 2 where kp and 
are the Fermi wave vector and energy, respectively, and L is the particle size. 
First, using an exact solution method for BCS Hamiltonian with finite num- 
ber of energy levels, for the well studied non-degenerate case we find a new 
feature. In that case, parity effect exhibits a minimum instead of a monotonic 
behavior. For d-fold degenerate states we find that the ratio of two successive 
parity effect parameters A p is nearly 1 + 1/d. Our numerical solutions for the 
exact ground state energy of negative U Hubbard model on a cubic cluster 
also give very similar results. Hence we conclude that parity effect is a general 
property of small Fermi systems with attractive interaction, and it is closely 
related to degeneracy of energy levels. 
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I. INTRODUCTION 



Recently, small metallic grains, with sizes down to 5-10 nm, became available for exper- 
imental study Superconducting transition in such grains manifests itself through the 
parity effect, namely dependence of electronic properties of the grain on whether the number 
of electrons in the grain is even or odd. Tinkham, Hergenrother and Lu Jl| showed that in 
a single-electron tunneling (SET) transistor [Q with a superconducting island, conductance 
channels open below a certain temperature T* < T c only for an odd number of electrons 
on the superconducting grain, whereas above T* they open symmetrically at odd and even 
number of electrons. This means that energy of the state with odd number of electrons 
shifts up by an amount A p || , so called parity effect parameter 

\ = E^-\{Ef N) +Ef N ^). (1) 

Several authors f||-|9| developed theories attempting to calculate properties of supercon- 
ducting state of the grain with discrete electronic levels. It was assumed that the level 
degeneracy is totally removed, and the only parameter distinguishing the small grain from 
the bulk sample is the ratio of the average level spacing S in a sample to the bulk supercon- 
ducting gap A, 5/ A. 

There are two main motivations for studying the effect of degeneracy. First, small su- 
perconducting particles can in principle be prepared in a perfect symmetric shape |L0] in 



which case the mesoscopic superconductivity parity effect will show novel features due to 
level degeneracy. The second and more important is the necessity of distinguishing between 
the level spacing due to size effect 8\ ~ hvp/L, and the average level spacing S 2 ~ ep/N 
which is much smaller than Si (N is total number of conduction electrons in metallic grain). 
The ratio §1/62 is of the order N 2 ^ 3 , which is very large for a typical mesoscopic particle of 
size 0.1/mi corresponding to N ~ 10 6 . The realistic energy structure of small grain is there- 
fore not equidistant spectrum with level spacing Si, adopted in most papers on mesoscopic 
superconductivity, but rather peaks in the density of states with characteristic distance Si 
accompanied by finer structure with smaller energy separation S2 ■ To understand which 
effects such finer structure may have on superconductivity in mesoscopic particles, we adopt 
a model of 5i-spaced levels with the degeneracy d (in the rest of our work we drop the index 



and denote level spacing by S). Calculation of superconducting condensation energy flT 
shows that energy binding per particle is not A, as may be expected from a naive picture of 
particle binding resulting in energy decrease 2A, but rather of order A 2 ftp. Adding extra 
two electrons to superconductor decreases energy by 2A but at the same time provides a 
(small) shift of the chemical potential in such a way that the net energy change per particle 
is only A x A/ep. A realistic calculation of the parity effect must treat system of interacting 
electrons self-consistently as it is done in the standard BCS theory of bulk superconductors 
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In the theory of superconductivity of Bardeen, Cooper and Schrieffer jL^j , the coupling 
Hamiltonian is introduced in the form 

Hint = g C kT C -k| c -k'jC k 'T (2) 
kk' 
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where operator c kcr creates an electron in a state with momentum k and spin projection 
a =|, I- k and — k states are selected in expense of all other states (such as k and q with q 7^ 
— k), because only these two time-reversed states lead to singularity which results in Cooper 
instability identified by logarithmic divergence of the scattering amplitude \yi{2Tioje>/ t) as 
e — > fLl] . Pairing of non-time reversed states can be treated perturbatively, or can be 



included as a renormalization in the Fermi liquid picture. 

Assume that the single electron states iprj in metal are decomposed as 

= E °k>k (3) 
k 

where k are plane waves. We are going to denote the time reversal state corresponding to 
if) v by iprj which is given by the complex conjugate of i[) v . In general ip^ and ipfj are different 
states. Let us introduce the operators 

4 = E«k4- (4) 

k 

These are Fermi operators since c kcr and Cko- obey the Fermi statistics. Let us consider the 
sum = J2 V c l^ c li> which can be represented as 

^ = E(E«)4A- ( 5 ) 

kk' V 

If we form an interaction Hamiltonian using A^ and A as Hi n t = gA^A, H contains a singular 
part which is identical to Eq. ^| because of the identities a k = (a^ k )* and J2r) l a k| 2 — 1 (f° r 
real ip^, \rja) and \fja) turn out to be the same state and if there is no further degeneracy 
due to other symmetries H int reduces to the "toy model" of superconductivity which has 
been studied extensively in the literature [0). Therefore, the interaction Hamiltonian in a 
grain can in general be written as 

H in t = ^EcWlWVt- ( 6 ) 

We will further split index 77 into j, a where j denotes the energy levels and a denotes the 
degenerate states for this level. The only symmetry expected to hold in an irregularly shaped 
ultrasmall grain is the time-reversal symmetry which may lead to double degeneracy for the 
energy levels, for example current carrying states. On the other hand, for certain energy 
spectra and more generally for lattice symmetries, energy levels are strongly degenerate near 
the Fermi energy. For very small particles, splitting due to weak disorder or asymmetry is 
much smaller than the splitting due to size effect. In this case the system can be modeled as 
almost degenerate with two dimensionless parameters governing superconducting transition, 
namely S/A, the ratio of level spacing to bulk gap, and d, the degeneracy of split levels. 



II. PARITY EFFECT FOR DEGENERATE LEVELS 

In our study, for simplicity, we assume all energy levels are <i-fold degenerate and levels 
are equally spaced with level spacing being equal to 5. We start with the following pair 
interaction Hamiltonian which is properly modified to take the degeneracy into account 
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xs 

H = S^JCjaaCjaa- ~J £ 4^4^'^' i C J'»'h ( 7 ) 
jaa j,j'£S;aa' 

where second sum is over the set of levels S which are lying within the 2tujr> shell centered 
at Fermi level, that is j,j' G S = {— n c , . . . , n c } with n c being equal to integer part of hcuo/S, 
and a — 1, ... ,d runs over all the degenerate states for a given energy level. Here, [j'af j.) 
and \j'a' |) are time reversed states for which the matrix elements are much larger then all 
others [14]] . For example, for a grain where eigenstates are labeled by crystal momentum k, 
the two states are |k f) and | — k |). Note that there is another similar but different pair 
formed by |k j) and | — k f). m usual BCS problem, since there is a summation over k, 
both pairs are properly taken into account. On the other hand, when we sum over energy 
levels rather than individual states we must be careful in including both pairs. However, 
the model without double degeneracy can still be used to describe the superconductivity in 



systems with real wave functions [|T3 



Next, we introduce the boson operators 
and 

bj a = CjalCjap (9) 

It is easy to verify that bj a and bj a satisfy "boson" commutation relations 

[bja, b), a ,} = SjfSaar^ - 2n ja ). (10) 

Presence of the number operator on the right hand side is direct consequence of the Pauli 
principle. With these new operators, the pairing Hamiltonian can be rewritten as 

H = 28Y,jn ja -^r E blb fa , (11) 

ja j,j'£S;aa' 

This Hamiltonian can be split into two parts as H = Hi + H e ff where 

H 1 = 25 E 3 n ja 

j?S;a 

H eff = 25 Y, 3 n i« tfabj'c*', (12) 

j£S;a j,j'GS;a,a' 

and the sum in Hi is over non- interacting particles while in H e ff both of the sums are over 
interacting particles. Interacting particles are those which occupy levels inside 2huo shell 
about the Fermi level except the singly occupied level. Since the interaction is only between 
pairs, this singly occupied level is removed form the set 5*. This is the so called "blocking 
effect". As Hi and H e jf commute, eigenstates of H are product of eigenstates of Hi and 
H e ff, and eigenvalues of H are sums of corresponding eigenvalues of H 1 and H e ff. H 1 simply 
represents non-interacting particles in a potential, thus determining the ground state energy 
of H reduces to solving for H e ff. 
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The effective Hamiltonian H e ff has been studied with exact diagonalization methods by 
Mastellone, Falci, and Fazio (MFF) ||, and by Berger and Halperin (BH) |J for the non- 
degenerate case and with relatively small number of levels taken into account (n c < 15). 
However, there exists an exact solution for H e ff due to Richardson and Sherman (RS) ||15|| . 
This long forgotten exact solution has been re-introduced to condensed matter community 
by Braun and von Delft Jl6| in the context of ultrasmall superconducting grains. Here, we 
briefly review the Richardson-Sherman solution. First, H e ff is treated without the Pauli 
exclusion principle, that is hardcore Bose particles are treated as normal Bose particles. 
After diagonalizing that Hamiltonian for normal bosons by a canonical transformation, the 
following boson wave function is obtained 

Uh...jN)=Y.nfi w . 1 F h as) 

v k=l Zd ^ - h Pk 

where N is the total number of bosons, and in J2 P > P means summing over N\ permutations 
of indices pi, . . . ,pn, and the corresponding energy is 

E = E P1 + ... + E m . (14) 

Next, Pauli exclusion principle is imposed with the condition that ip — if any two ji are 
not distinct. This restriction is introduced by writing the wave function as 

ip(ji ■ ■ ■ Jn) = 9(ji ■ ■ ■ 3N)(j>{h ■ ■ ■ 3n) (15) 

where 9(ji . . .j N ) is equal to 1 if all jj's are distinct, and is equal to otherwise. Then, after 
some calculation an effective Schrodinger equation for is obtained 



X5 N N 

(25 ji + ... + 25 j N - E)4>(jt . . .j N ) - — ^2 XJt 1 - J2 S ik)(f>{ji ■ . ..Jn) = 0. 

(16) 

For a multiply degenerate single level, ground state energy is given by 

E = 2Ne- —N(d- N + 1), (17) 

CI 

where d is the pair degeneracy of the level. Next we focus on degenerate levels in the 2huj£> 
interacting shell leading to a set of coupled equations for iV pairs and M levels which are 
d-fold degenerate: 

d p^g ^p-^q d a j=-n c 2 °3 ~ &q 

where 

E p ^ E q , for p q. (19) 
The ground state wave function is given by ([13]) up to a normalization constant and the 



corresponding ground state energy is calculated from (|14|) . 
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At some values of X5/d, depending upon the state under consideration, Eq. [18] has 
singularities. At a singularity the restriction (|19|) is not satisfied. However, the domain of 
validity of the solution can be extended by letting E q to be complex. Complex roots E q 
occur in complex conjugate pairs 

-S27 = x-y + iy-y, -E27— 1 = ^7 — W^i 

7 = l,...,iV/2, (20) 

where x 1 is real and y 7 can be either real or pure imaginary. It turns out that if N is not 
even, one of the roots remains real for all A. This form of E q preserves the reality of the 
ground state energy E (Eq. [L4D, and also the reality of the ground state wave function |L5 



At a singular point, no more than two roots can be equal, and these two equal pair energies 
are both 25 ji for some ji. The desired roots of Eq. [18] for the ground state should satisfy 



lim E a 



q-1 

d 



1)5, g = l,..., N, (21) 



that is, in the non-interacting system the lowest N levels are occupied by N pairs while 
levels from N + 1 to M are unoccupied ([...] denotes integer part of the division). We solve 
Eq. [TS| by a globally convergent Newton-Raphson method |T7| for complex E q ( |2"0"D with the 
conditions implied by Eq. [2~I| . 

For the non-degenerate case singularities of Eq. [18] can be removed by a suitable change 



of variables ||18[| . However it is not possible to generalize that method to degenerate case. 
In this case roots are complex for all values of A (if N is odd, one root remains real for all 
A). Fig. [I] shows typical behavior of roots for the double degenerate case (see Richardson 
|18| for the behavior of roots in the non-degenerate case). 

The ground state energy E^ for a given number of electrons n allows us to calculate 
the parity effect parameters 

_ |^n(2n c d+i-l) _ ~^_^(2n c d+i-2) _|_ j^(2n c d+i) ^ | (22) 

where i = 1, . . . , 2d. Here by n we mean the number of electrons in the thin shell around 
the Fermi level participating the pairing interaction. We assume that the shell is composed 
of n c levels below and above the Fermi energy. Note that with this definition Eq can take 
2d different values and that is why i index running from 1 to 2d exhausts all possible values 
of the right hand side. 

In principle, to study the effect of finite energy level spacing, we should fix the Debye 
frequency u>d, which is assumed to be less affected by the boundary conditions, and change 
the number of levels n c . We estimate that typical n c values lie in the range of Tiujd/5 ~ 
50 — 2000, and for 5 ~ A, we have n c ~ 100. Numerical solution of Eq. [18| becomes 
more complicated with increasing n c values for degenerate cases. For this reason we use an 
alternative approach where we vary the Debye frequency up and the coupling strength A in 
order to vary 5/ A ratio. The motivation from physical point of view is the expectation that 
number parity effect, if it exists, should mainly be a function of dimensionless parameter 
5/ A. Our numerical results support our alternative approach, as well. Given a 5/ A ratio, if 
we repeat our calculations for increasing values of n c , we observe that Aw/A ratios do not 
change very much. 
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Exact diagonalization approach of MFF || and BH || for the non-degenerate case showed 
that both ground state properties and excitation gap are parity dependent and functions of 
the ratio of level spacing to the BCS gap, 5/A. However systems addressed by MFF and BH 
are limited to relatively small number of levels, n c , taken around Fermi level. Practically 
n c ~ 15 is an upper limit for any exact diagonalization scheme (either Lancsoz method or 
other methods) because of large memory space requirements. In Fig. ^| we reproduce the 
dependence of parity gap parameter, A«, upon level spacing 5. Here our aim in reproducing 
these results, which are obtained by using much larger values of n c (such as n c = 500) for the 
non-degenerate case, is to compare them with those of MFF and BH. When n c increases we 
obtain a different behavior of which was not observed by MFF and BH. As the number 
of levels n c increases, instead of a monotonic behavior, / A curve exhibits a minimum 
at 5/A ~ 0.5, which can not be observed with smaller n c values (see Fig. @(a)). On the 
other hand, for the dependence of upon 5 we observe the same behavior as MFF and 
BH, that is a minimum at 5/A ~ 1 (see Fig. 0(b)). Von Delft and Braun [|19j show similar 



results for A^ in the non-degenerate case in order to compare RS exact solution to earlier 
approaches to the problem of mesoscopic superconductivity. 

Next we performed calculations for degenerate cases (d > 2). In Fig. |3| we present the 
results for doubly degenerate case (d = 2) which show that parity effect is still there. Similar 
results are obtained for higher degeneracies (d > 2). Some important conclusions can be 
drawn by comparing Fig. ^] and Fig. [|, immediately. Firstly, A^s repeat themselves with 
a periodicity of two for the non-degenerate case, whereas they have a periodicity of four for 
the doubly degenerate case, that is the parity effect parameters repeat themselves with a 
period of 2d. Secondly, one can immediately observe that there is a remarkable difference 
between the ratios of two different A^s for the degenerate and non-degenerate case. For 
example, the ratio A^/A^ (Fig. |(a) and (b)) at a fixed 5/A value (~ 10), is about 5, 
while A^ / Aip for doubly degenerate case (Fig. ||(b) and (c)) is about 1.5 for the same 
value of 5/A. We will come to a rough estimation of these values in a short while. 

As it has been mentioned earlier, for certain energy spectra and more generally for 
lattice symmetries, energy levels are degenerate near the Fermi energy. In case of a d-fold 
degenerate single level with our notation n c = 0, and the ground state energy ([H]) measured 
from the Fermi level is given by 

E^ = ^[N*-(d + 1 + { - ir )N] (23) 

where N is the integer part of n/2 as above. Parity effect appears in the second term on the 
right hand side. Depending upon whether n is odd or even, the factor in front of N becomes 
d or d + 1, respectively. This is nothing but the blocking effect of the single electron. It is 
clear that when d becomes larger, this effect will become less important. By substituting 
Eq. [23| into Eq. [22| we obtain 

A (odd) i 

+ t (24) 



A (even) d 



where odd (even) stands for i odd (even) in Eq. Ell Note that we consider a fixed chemical 



potential /z in derivation of (p4[) . Therefore this ratio (Ej|) is valid for a parity effect parameter 
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for which all ground state energies, Eq , used in Eq. |2| calculated with the same 
chemical potential /x. Comparing Figs. 0(b) and (c) at 5/ A ~ 1 we find that A^'/A^' ~ 
1.33, which is quite close to 1 + 1/2 = 1.50. Moreover, for 5/ A ~ 10, we find A^/A^ ~ 1.51 
which is even closer to the ratio above. For 4-fold degeneracy, where 1 + 1/d — 1.25, we 
obtain similar ratios. For example A^/A^ ~ 1.17 at 5/ A ~ 1, and A^/A^ ~ 1.24 
at 5/ A ~ 10. These results are not unexpected since as 5/ A increases, contribution of 
electrons at the Fermi level to the ground state energy becomes dominant which lets single 
level result become a better approximation. On the other hand when /i shifts, Eq. [3] is 
not valid anymore and the ratio A^ M) /A^ ven) becomes very different from 1. For example, 
A^/A^ ~ 5 for non-degenerate case and A^/A^ ~ 5 for doubly degenerate case at 
8/ A ~ 10. 



III. PARITY EFFECT IN ATOMIC CLUSTER 

We supplement the investigation of a parity effect "from above" (from macroscopic sizes 
down to mesoscopic sizes) by an investigation "from below", that is, starting from small 



clusters of atoms coupled at sites by some interaction energy ||20|| . Unlike Eq. y, we choose 
the negative-?/ Hubbard interaction 



N 

H = -t^aL a ja + Uj2 n ^ n H (25) 

where i and j numerate atomic sites, a\ a is the creation operator at site i with spin projection 
a, and 7ii a is the number operator. This above Hamiltonian will reduce to a BCS Hamiltonian 
(§) in the limit of large number of sites and small U with the difference that the interaction 
( P5p is a non-retarding one, and the interaction shell (the cut-off energy) at Fermi energy is 
increased up to the value of the order of Fermi energy itself. In the absence of interactions for 
energy levels of cubic cluster, we find — 3t, —t, t and 3t where —t and t are triply degenerate, 
and hence level spacing S is given by 2t. An exact solution of Hamiltonian (1251) in a cubic 



cluster allowing for maximum of 16 electrons has been studied earlier ||21|| . The Hamiltonian 
matrix corresponding to Eq. [25| has a maximum dimension of ~ 6 x 10 4 after reduction of 
the size with symmetry consideration, and ground state eigenvalue is calculated by an exact 
diagonalization method (of non-Lanczos type) ||22|| . Fermi operators aj CT are represented as 



a 2 | = u®u®a<8)'U<g>M...<g>M (26) 



a 8 i = f®f®f®t>(g>t>...Cg>a 
where a, u, v are 2x2 matrices 






(27) 



S 



and ® is the Kronecker matrix product symbol. By changing U, we change the bulk energy 
gap A which we calculate by using BCS gap equation. For —1 > U/t > — 10, approx- 
imating the density of states by 1/12£, we see that the coupling parameter A ~ U/Ylt 
changes between 0.08 and 0.83. Therefore, at least for small \U\ values, we are in the weak- 
coupling BCS regime where approximate solution of the gap equation can be written as 
A = 12i exp(— 12t/\U\). Non-zero parity effect parameter, A p , is directly seen from the 
energy versus electron number plot (see Fig. 4). We observe that Eq is not a monotonic 
function of n. The ground state energy exhibits dips for even number of electrons. This is 
an unambiguous indication of the parity effect, which is a direct manifestation of Cooper 
pairing. 

The inset in Fig. 4 shows that = E^ — (E^ + Eq 8 * > )/2 is non-zero as long as 
there is an attractive interaction at sites, i.e. U < 0. We find, on contrary, that A p = at 
U > and hence conclude that repulsive Hubbard model does not lead to superconductivity. 
Nevertheless in principle this result does not exclude the possibility of superconductivity in 
a larger system with U > 0. 

Presence of parity effect for negative-?/ Hubbard Hamiltonian is important since it shows 
that dependence of superconducting properties on whether there are even or odd number 
of electrons is not specific to BCS Hamiltonian. Although our rule that the ratio of two 
successive parity effect parameters is 1 + 1/d does not work, we clearly see that A W values 
form groups in parallel to the degeneracy of levels. For example, the first line in Table | 
corresponds to a jump from —3 state to —1 state, while the next five lines correspond to 
energy levels —1. It is possible to understand different behavior of A^ for % = 3, 9 and 
15 by introducing analogue of chemical potential given by // = (-Ehomo + -E'lumo)/2 where 
HOMO and LUMO stand for highest occupied and lowest unoccupied molecular orbital, 

(k) 

respectively. For i = 3, 9 and 15, /i is different in each Eq value used in calculations of 
A». For example, in case of = \E^ S) - (E^ + E^)/2\, /i = 0,-1 and 1 while for 
A^ 11 ) = I-Eq 10 ^ — + Eq )/2|, fi = 1 for all three ground states. Figure [| shows three 
different parity effect parameters, A^\ A p 10 ^ and A^ 11 ) as a function of level spacing 5 = 2t. 
It is remarkable that curves are very similar to those obtained for BCS Hamiltonian. Again, 
due to degeneracy, A^ 9 ) which involves the jump — 1 — ► 1 exhibit a slightly different behavior 
in comparison to A^, 10 - 1 and A^ 11 ). As far as filling of molecular orbitals is concerned i = 9 
case is analog of Fig. 0(a), while i = 10 and 11 correspond to Fig. |3|(b) and (c), respectively. 

IV. CONCLUSION 

In conclusion, we have studied the effect of degeneracy of discrete energy levels on su- 
perconducting properties of a small metallic grain. We observe that parity effect parameter, 
which is a measure of the dependence of energy on whether the number of electrons in the 
grain is even or odd, exhibits a behavior similar to the non-degenerate case for small d. 
We reproduced the behavior of parity effect parameter in the well studied non-degenerate 
case in order to compare our exact results with previous work. In that case the parity ef- 
fect parameter exhibits a minimum instead of a monotonic behavior. For d-fold degenerate 
states it turns out that there are 2d different parity gaps and furthermore both approximate 
analytic solutions and exact numerical results suggest that ratio of two successive parity 
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effect parameters is nearly 1 + 1/d. Therefore careful measurements of parity parameters 
can be used to determine the degeneracy of energy levels. With increasing degeneracy, the 
pairing effect is suppressed. Although there is no direct evidence for existence of SET tran- 
sistors with a metallic grain of perfect crystal structure, observation of perfect geometries 
in small grains prepared by vapor condensation technique opens a possibility for having 
highly symmetric structures. As we discussed above smallness of the average level spacing 
in comparison to the level spacing due to finite size effect leads to a nearly degenerate energy 
spectrum. Convergence of the ratio of two successive parity effect parameters to 1 + 1/d as 
the level spacing increases, can be used to detect the degeneracy. 

We also show that parity effect is not specific to BCS Hamiltonian by solving negative- 
U Hubbard model for a small atomic cluster exactly. We can say that this is the most 
rigorous way to treat the problem since it does not involve any approximation like pairing 
assumption. Our results clearly show that parity effect exists and hence we conclude that 
it is model independent and a general property of small Fermi systems with attractive 
interaction. Furthermore, A® /A curves exhibit a behavior which is very similar to BCS 
case. Finally, grouping of AW values according to energy levels of atomic cluster shows that 
degeneracy still manifests itself. 

This work was partially supported by the Scientific and Technical Research Council of 
Turkey (TUBITAK) under grant No. TBAG 1736 
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TABLES 



TABLE I. A 



(0 



E, 



(i-i) 



-2) 



for different values of i with U 



-1. We 



observe three relatively large A p ^ values (i = 3, 9, 15) corresponding to jumps between the energy 
levels of the cubic cluster (—3 — > —1, —1 — > 1, and 1 — >• 3). 







3 


1.003 


/i 

4 


n on/i 
U.ZU4 


5 


0.212 


6 


0.201 


7 


0.206 


8 


0.197 


9 


1.010 


10 


0.197 


11 


0.206 


12 


0.201 


13 


0.212 


14 


0.204 


15 


1.003 
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FIGURES 



FIG. 1. Behavior of roots of Richardson-Sherman formula in the doubly degenerate case, (a) 
For even number of bosons, N, all roots are complex conjugate pairs, and real parts, x 7 , are plotted 
with respect to A. Since the roots, £27 and -E27-1 are complex conjugate pairs, their imaginary 
parts (y 7 ) cancel each other and the ground state energy becomes twice the sum of all x 7 . (b) 
When the number of bosons, N, is odd, N — 1 roots E q are complex conjugate pairs and one 
remaining root (E21, shown by a darker solid line in (b)) is real for all A. Again only the real parts 
of roots, x 7 , are plotted with respect A. 



FIG. 2. Parity effect parameters as a function of level spacing for non-degenerate case, (a) and 
(b) n c = 15, 30, 60, 120, 240, 360, 500 from bottom to top. (a) With larger values of n c , which are 
more realistic in comparison to experiments, we obtain a behavior of A^ which was not observed 
previously. Instead of monotonically decreasing towards Ap^/A = 1, the curve makes a minimum 
at 5/ A ~ 0.5. The inset shows the details around the minima. However it still approaches to 
one after this minimum as expected, (b) A p 2 ^ has the same expected behavior as shown by many 
previous studies. It has a minimum at 5/ A ~ 1 and after this minimum it turns upward towards 
A p 2) /A = 1 



FIG. 3. Parity effect parameters as a function of level spacing for doubly degenerate energy 
levels. n c = 20,25,30 from bottom to top in (a), (b) and (d), from top to bottom in (c). The 
immediate fact observed by comparing this figure with the one for the non-degenerate case exhibits 
itself when the ratios of different A^s are considered. For example, when we take ratio of A^ 
to Ap 3 - 1 (compare (b) and (c))we see that its value is much smaller than such a ratio in the 
non-degenerate case (compare Fig. 2(a) and (b)). 



FIG. 4. Dependence of ground state energy Eq (in units of t) upon number of particles n for 
U/t = —2, —3, —4 from top to bottom. The ground state energy exhibits drops for even number 
of electrons. Lower inset shows the cubic cluster on which negative-?/ Hubbard model is solved. 
Parity effect parameter A p versus on-site interaction U at half filling is shown in the upper inset. 



FIG. 5. Parity effect parameters A p 9 \ A p 10 ^ and Ap 11 ^ as a function of level spacing 5 = 2t, for 
cubic cluster where (a) A p 9) = \E^ 8) - (E^ + E^)/2\ (b) A p 10) = \E^ 9) - (E^ + ^ 10) )/2| and 
A P U) = \E^ 0) - (£j 9) + E^ U) )/2\. The insets show the corresponding energy levels of the cluster. 
Since A^ involves the jump — 1 — ► 1, it is relatively large in comparison to Ap 10 -* and Ap 11 -* as it 
can be seen in Table I. 
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Fig. 4 
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Fig. 5 
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